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I A systematic method is presented for constructing effective Hamiltonians 

\^ ' for general phonon-related structural transitions. The key feature is the ap- 

o' 

^D ' plication of group theoretical methods to identify the subspace in which the 






effective Hamiltonian acts and construct for it localized basis vectors, which 
are the analogue of electronic Wannier functions. The results of the sym- 
metry analysis for the perovskite, rocksalt, fluorite and A15 structures and 
the forms of effective Hamiltonians for the ferroelectric transition in PbTiO^ 
and BaTiO^, the oxygen-octahedron rotation transition in SrTiO^, the Jahn- 
\—i ' Teller instability in Lai-x{Ca, Sr, Ba)xMnOs and the antiferroelectric transi- 

c^ ! 

tion in PbZrO^ are discussed. For the oxygen-octahedron rotation transition 
in SrTiO^, this method provides an alternative to the rotational variable ap- 
proach which is well behaved throughout the Brillouin zone. The parameters 
appearing in the Wannier basis vectors and in the effective Hamiltonian, given 
by the corresponding invariant energy expansion, can be obtained for individ- 
ual materials using first-principles density-functional-theory total energy and 
linear response techniques, or any technique that can reliably calculate force 
constants and distortion energies. A practical approach to the determination 
of these parameters is presented and the application to ferroelectric PhTiO^ 
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discussed. 
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I. INTRODUCTION 

Interest in the construction of effective Hamiltonians to relate macroscopic properties of 
materials to microscopic information has recently increased dramatically. Construction of 
an effective lattice Hamiltonian is greatly facilitated by the use of a localized basis for the 
description of crystal structure distortions. For example, structural phase transitions in per- 
ovskite and rocksalt structure systems have been studied using model Hamiltonians with a 
"local mode" basis to describe the unstable lattice modes relevant to the transition. While in 
earlier work the microscopic model parameters had to be determined empirically,aa recent 
advances in the accuracy and efficiency of first-principles bandstructure and total energy 
methods have made it possible to calculate these parameters directly from first principles, 
and encouraging agreement with experimentally observed transition properties has been 
obtained.Q In addition, the study of electronic properties of systems which exhibit lattice 
instabilities coupled to the electronic states is facilitated by the description of the lattice 
in a localized basis, putting it on an equal footing with the usual tight-binding description 
of the electronic states.^ However, the approximate "local mode" basis is expected to be 
inadequate for high-accuracy first-principles investigations. Furthermore, explicit construc- 
tion of a localized basis set and optimal effective Hamiltonian has proved to be especially 
problematic for complicated crystal structures and structural transitions involving modes at 
the zone boundary or interior of the zone. 

In this paper, we present a systematic procedure for the construction of an effective lat- 
tice Hamiltonian in a localized basis set, applicable to arbitrary systems and structures and 
specifically designed for practical use in first-principles studies. This approach is suitable 
for the investigation of a general phonon-related structural transition^ involving a small 
change in structure which can be generated by freezing phonon displacements into a high- 
symmetry prototype structure. We implement this procedure for structural transitions in 
the perovskite, rocksalt, fiuorite and A15 structures. Group theoretical methods, which since 
Landauty have been central to the phenomenological study of structural phase transitions. 



are essential to our development of a procedure which is systematic and easy to implement 
for arbitrary structures and transitions. The resulting localized basis vectors are the lattice 
analogueQBO'y of electronic Wannier functions.E^'EJ These Wannier basis vectors span an 
invariant subspace containing the low-frequency normal modes relevant to the transition. 
This subspace is energetically decoupled at quadratic order from the higher frequency eigen- 
modes, permitting the construction of an effective Hamiltonian for the transition including 
only lattice distortions in the subspace. The model parameters for a particular material 
can be obtained using any method that can reliably calculate force constants and distor- 
tion energies. For the force constant calculation, the most efficient high-accuracy method is 
density-functional-theory (DFT) linear response .cJlli As we will see, once the force- const ant 
matrices have been calculated, an effective Hamiltonian can be constructed with very little 
additional effort. The final result is a model system which reproduces the finite-temperature 
transition behavior of the original lattice Hamiltonian, while the simpler form and reduction 
in the number of degrees of freedom per unit cell make it suitable for study by methods such 
as mean field theory and Monte Carlo simulation. 

In the next section, we present the derivation of the effective Hamiltonian in a sym- 
metrized localized (Wannier) basis. In Section HI, we discuss in detail the construction 
of the Wannier basis and resulting effective Hamiltonian for the simple example of the di- 
atomic chain. In Section IV, this construction is generalized to arbitrary structures and 
applied to the structural transitions in perovskite compounds. In addition, the results of 
the construction for the rocksalt, fiuorite and A15 structures are briefiy described. Section 

V is a discussion of the practical implementation of this approach in a first-principles frame- 
work, illustrated by application to the ferroelectric transition in PbTiOs. Finally, Section 

VI summarizes and concludes the paper. 



II. FORMAL DERIVATION OF THE EFFECTIVE HAMILTONIAN 

Consider a crystal with Uat atoms per unit cell and N unit cells with Born-von Kar- 
man boundary conditions. Each possible configuration of the crystal is fully described by 
specifying the generalized displacement of each ion (i.e. the displacement multiphed by the 
square root of the ion mass) relative to a fixed high-symmetry reference structure, which 
should be chosen so that the physically important configurations are obtained by small dis- 
tortions. In general, this will be the same as the structure of the high-symmetry phase. 
The set of all ionic displacement patterns forms a Sn^fA^- dimensional vector space, called 
the "ionic-displacement space." In the standard basis, in which each basis vector is the 
unit displacement of a single ion in a single cartesian direction, the coordinates of a par- 
ticular configuration are the generalized displacements {uja, j = l.UatN; a = x,y,z}. 
Another basis for the ionic-displacement space is obtained from the normal modes {eg^; k G 
Brillouin zone, A = 1, Suat}- In this basis, the coordinates {cj:^} are complex numbers, with 
the understanding that the ionic-displacement pattern is obtained by taking the real part 
of the vector J2k xCj^x^kx- '^^^ (^a-a) transform according to irreducible representations of 
the space group of the high-symmetry reference structure Qo- For example, the action of 
translation through a lattice vector is (9^ eg^ = exp{ik ■ R)e^y 

For the construction of an effective Hamiltonian, we perform a decomposition of the 
ionic-displacement space into invariant subspaces, called "band subspaces," spanned by one 
or more entire branches of the normal modes and closed under QS^ Each band subspace in 
this decomposition should be minimal in the sense that it cannot be written as the direct sum 
of two band subspaces. The number of band subspaces is Ug and they are indexed by A, with 
the number of branches in each subspace nA, indexed by / (note that I]A=i,ns "-a = 3r2aj). 
n\ will also be referred to as the multiplicity of the band subspace. 

For each band subspace A, we seek a localized basis set with definite symmetry properties 
with respect to Qq. These properties are established by choosing a Wyckoff symbol for Qq, of 
multiplicity n^, and an /^-dimensional irreducible representation (irrep) of the corresponding 



site symmetry group.lljOEl The choice of Wyckoff symbol and site symmetry group irrep is 
constrained by the requirement that the representation of Qq defined by the locahzed basis 
have the same decomposition into irreps of Qo as the particular band subspace A for which a 
basis is being sought. Specifically, a set of basis vectors {u^Aii_R„i ^ ~ -'-' ^^l ^^ associated with 
a particular site in the crystal (the ith Wyckoff position added to the lattice vector Rn) and is 
a set of partner functions of the chosen irrep of the site symmetry group. These basis vectors 
are localized, with ionic displacements decaying rapidly with distance from the associated 
site. The basis vectors associated with the other sites of the chosen type are generated by 
acting on the given basis vectors with symmetry operations from Qq. For example, subsets 
of localized basis vectors are related to each other by translation through a lattice vector 
'^MjRn+R ~ ^R '^AijRn' Clearly, the product of the number of Wyckoff positions n^ and 
the dimension of the site symmetry group irrep Ig must equal the multiplicity of the band 
subspace n^, and we replace the former two indices by a combined index / = l,nA. The 
decomposition of the locahzed basis representation of Qq into irreps of Qo can be computed 
from formulas such as Equation 14 in Ref. O or by computing, for various k, the character 



table of Qj^ in the subspace of Bloch sums. Usually, the appropriate choice of Wyckoff symbol 
and site symmetry group irrep can be determined just from the space group irrep labels of 
the eigenmodes at the high-symmetry points. 

The basis vectors from each band subspace in the decomposition can be combined to form 
a basis for the full ionic displacement space. In this basis, the coordinates of a particular 
configuration are the real numbers {^^jj^.; A = l,ns; I = l,n/^; i = 1, A^}, which determine 
the ionic-displacement pattern as the vector J^ah^air'^air ■ These localized basis vectors 
are the lattice analogueu'BOa of electronic Wannier functions.ta'EJ The symmetry properties 
of the Wannier basis vectors lead to a convenient form for the Taylor expansion of the 
lattice energy in the corresponding coordinates, allowing easy identification of independent 
parameters. Furthermore, the expansion has the important feature that it contains no 
quadratic cross terms between different band subspaces. Extra considerations apply when 
the various branches of normal modes include some with polar character and/or are not well 
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separated in energy; these will be discussed in Section V. 

With this basis set, the derivation of an effective Hamiltonian proceeds as follows. The 
classical partition function for the lattice is: 



Z (X {duja}exp{-f3Hiat{{uja})) (2.1) 



where j runs over the HatN atoms, the electron degrees of freedom have been eliminated 
through the Born-Oppenheimer approximation, and the generalized ionic displacements 
{uja} are measured relative to a high symmetry reference structure, as described above. 
The momentum integral leads to a trivial w-independent factor and is not explicitly written 
here. 

It is straightforward to change variables in Tiiat from standard basis coordinates {uja} to 
Wannier basis coordinates {'^Ai^R.} ^-i^d Taylor expand Tiiat around the reference structure. 

(2) 

By construction, the quadratic part Ti^J decouples into a sum over the band subspaces: 



A=l,ns 

where 



'^A ii^AlRj) ~ A^ A^ '^lil'j^AlRi^Al'Rj- 

l,l'=l,nA i,j=l,N 

The next step is to choose a band subspace or direct sum of band subspaces as the 
symmetry-invariant subspace in which the effective Hamiltonian will act. This subspace 
should be as small as possible, yet should include all modes which have higher-order terms in 
the expansion oi Hiat which are important in describing low-energy configurations, and which 
therefore cannot to a good approximation be dropped from the expansion in the evaluation 
of the partition function near the transition. For the calculation of thermodynamic behavior 
at temperatures in the vicinity of the structural transition, it is essential to include higher 
order terms for the branches containing the unstable modes associated with the transition. 
Strain degrees of freedom may not have important anharmonicity in themselves, but still 
have physically important higher-order coupling to the unstable modes, and thus should be 
considered for inclusion in the subspace. Here and in the rest of this section, we write the 
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expressions for the case that the effective Hamiltonian acts in a single band subspace Aq, 
since generahzation to the case of the direct sum Ai © A2 © ... © A^ is straightforward. The 
correct choice of subspace permits the following approximation to Hiat- 

-Hlat ^ K//({eA„,^J) + E '>it\{iMR.}) (2.2) 

A^Ao 

Within this approximation, anharmonic terms appear only in Tie// and, by construction, 
there are no cross terms between Aq and A ^ Aq. This has the important consequence 
that the integration over the degrees of freedom not included in the effective Hamiltonian 
subspace can be performed analytically. Specifically, the partition function integral factor- 
izes into an integral over Aq and another over A 7^ Aq. The integral over A 7^ Aq is an 
elementary Gaussian integral which leads to a temperature-dependent factor unimportant 
for the structural transition. The remaining integral has the form of a partition function 
for a system of nAo-component spins on a lattice {-Rj}, which recovers the original partition 
function Z: 

j{dU,mM^Pi-meff{{U,iR.})) ^ Z. (2.3) 

An explicit expression for 'Heff is obtained by making a Taylor expansion in symmetry 
invariant combinations of the {^^ ^^ }. This expansion is greatly facilitated by the symmetry 
properties of the Wannier basis. An expression with a finite number of undetermined param- 
eters is obtained by setting expansion coefficients beyond a certain range or order to zero or 
writing them in a parametrized form. The resulting spin-like model is suitable for statistical 
mechanical analysis via mean field theory, Monte Carlo simulation, or other appropriate 
method. The values of the expansion coefficients are different for different materials and are 
to be obtained from first principles calculations, as will be discussed in Section V. 

III. EFFECTIVE HAMILTONIAN CONSTRUCTION: THE DIATOMIC CHAIN 

In this section, the explicit construction of the effective Hamiltonian in a localized Wan- 
nier basis will be presented for the simple example of the diatomic chain.EJ This calcula- 



tion demonstrates most, though not quite all, of the important principles, including the 
determination of the symmetry properties of the Wannier basis vectors. The remaining 
generalizations for application to arbitrary structures will be made in Section IV. 

The high-symmetry reference structure of the "ionic" diatomic chain model we consider 
has lattice constant a = 1 and ions A and B at positions n and n + ^ {n = 1,..., A^), 
respectively. The potential energy in our model is given to quadratic order in {unAjUns}, 
the displacements from the reference structure, by 

^[7^{UnB - UnAf + ^^(Mn-l.B - UnAf + T^i'^n+l^B " Unsf] (3.1) 

n ^ ^ ^ 

as shown in Fig. |l]. In the rest of this section, we set the masses of the ions m^ = thb = 1 
and the ratio - = 0.25. 

a 

Analytical diagonalization of the dynamical matrix at each k in the first Brillouin zone 
yields two branches (A = 1, 2) with eigenvalues e!^{k) = uj\{k), shown in Fig. ^ The normal 
modes are labelled by the irreps of the space group {{U\nx}, U = E,!; n = 1,...,A^}. 
These are Tg and T^ at k = 0, Xg and X^ at A; = vr, and k otherwise. It can be seen that 
the ionic displacement space of the diatomic chain decomposes into two band subspaces 
A = 1,2, each spanned by the nondegenerate eigenmode branch of corresponding A. Since 
the number of branches in each subspace is one, we drop the index /. 

As we have seen in Section II, construction of an effective Hamiltonian involves the choice 
of an invariant subspace. This subspace should be a band subspace (or direct sum of band 
subspaces) which includes all modes relevant to the particular transition being modelled. In 
the diatomic chain example, we imagine that the relevant mode is X„. The band subspace 
A = 1 has the symmetry label Xu, while the band subspace A = 2 has the symmetry 
label Xg. Therefore, the subspace in which the effective Hamiltonian will act is A = 1. If 
instead the relevant mode were Xg, the corresponding effective Hamiltonian subspace would 
be A = 2. 

For the chosen effective Hamiltonian subspace A = 1 or A = 2, we seek a localized 
basis set labelled by a unit-multiplicity Wyckoff symbol and a one-dimensional irrep of the 



corresponding site symmetry group, since both band subspaces are of unit multiplicity. The 
space group of the diatomic chain has two Wyckoff symbols of unit multiplicity: la, (at 
x=0), and lb (at x=l/2). For both la and 16, the site symmetry group is 1, with two 
one-dimensional irreps Ag (even) and Au (odd). The space group representations defined by 
each of the four distinct labels {la,Ag), {la,Au), {lb,Ag) and {lb,Au) decompose at the high- 
symmetry k-points A; = and k = it into the space group irreps (r+,X_|_), (r_,X_),(r+,X_), 
and (r_,X-(.), respectively. Since the space group irrep labels of the localized basis set must 
match those of the band subspace throughout the Brillouin zone, we find that for both band 
subspaces the label of the Wannier basis is uniquely determined. For A = 1, it is {la,Au), 
and for A = 2, it is (16, A„). 

Since individual ionic displacements are a localized basis with the required symmetry 
properties for the full ionic displacement space, it might seem that the Wannier basis vec- 
tors would always be obtained from occupied Wyckoff positions combined with vector rep- 
resentations of the site symmetry group, as in the diatomic chain example. In fact, this is 
not the case. The simplest example in which this can be seen is the "molecular" diatomic 
chain.liil For the present discussion, we consider the closely related triatomic chain model in 
Fig. ^ The high-symmetry reference structure has lattice constant a = 1 and ions A at 
positions n and ions B at positions n ± | (n = 1, ..., N). The space group is the same as 
that of the "ionic" diatomic chain discussed above. The space group irreps for the normal 
modes at the high-symmetry k-points are F^, F„, and r„ at A; = and Xg, X„ and Xu 
at /c = TT. All space group irreps at the two points are compatible, so depending on the 
details of the potential, four types of branches defining unit-multiphcity band subspaces can 
appear. The Wannier basis of the band subspace (F„, X„) is labelled by (la,y4„), so that the 
basis vectors transform like ionic displacements at occupied Wyckoff positions. However, 
the Wannier basis of the band subspace (F^, Xg) is labelled by (la, Ag) and the Wannier 
bases of the two "mixed-label" band subspaces (F„, Xg) and (F^, Xu) are labelled by the 
unoccupied Wyckoff position 16, midway between the two B atoms, and the site symmetry 
group irreps A^ and Ag, respectively. For the triatomic chain, with an appropriate choice 
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of a and 7, it may happen that there will be two band subspaces of the (F^j, X^) type. 
Both will be described by Wannier basis vectors with labels (la,y4„), but the Wannier basis 
vectors for the two band subspaces will be orthogonal to each other. 

For the complete construction of an effective Hamiltonian, it is not enough to know the 
symmetry properties of the basis vectors; rather, an explicit expression for the Wannier basis 
vectors in terms of ionic displacements is needed. The basis vectors are determined by the 
potential energy and resulting set of normal modes which define the band subspace. If the 
normal modes are known at all fc, the basis vectors can be obtained exactly. For each A, the 
expression for wa^, the Wannier basis vector centered in unit cell n, in terms of {cfcA} is 

WAn = — Y^ exp{i(j)Aik))ekAexp{-ikn). (3.2) 

— 7r</c<7r 

An inverse transform can be performed to express the {cfcA} in terms of the Wannier basis 
vectors: 

CfcA = Gxp{-i(t)A{k)) Y^ exp{ikn)wAn- (3.3) 

n 

In evaluating Equation \i.2\, the phase function 0a (^) can be freely chosen to maximize the 



localized character oiwAn, with the only constraint being the requirement that the resulting 
Wannier basis vectors have the desired symmetry properties. For the diatomic chain, the 
normal modes are easily obtained at all k. We follow the convention that the real part of 
the 2-component eigenvector (A^^, A^^) of the 2x2 dynamical matrix at each k gives the 
displacements of the A and B ions in the unit cell n = 0. The displacements of the A and 
B ions in the unit cell n then are Re{A}^^exp{ikn)) and Re{Al.^exp{ikn)). A^^ and A^J 
can be chosen to satisfy |A^^p + |A^^p = 1, A^^* = A_^^ and A^J* = A_^^ and to be 
smooth functions of k, ensuring that the corresponding normal mode e^A smoothly depends 
on k. For A = 1, we choose (A,t^, Ajt^) such that Aj^^ is real and positive for all k, and 
4'A{k) = 0. The resulting Wannier basis vector win is shown in Fig. ^. 

Except in special cases such as the diatomic chain, knowledge of the normal modes at 



all k is not available and the expression in Equation |3.2| cannot be used to construct the 
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Wannier basis vectors. However, a good approximation can be obtained from knowledge of 
the eigenvectors at a subset of k-points through a fitting procedure which will prove to be 
practical for use with first-principles calculations. For a particular choice of Wyckoff symbol 
and site symmetry irrep, a parametrization of the displacement pattern of the Wannier basis 
vectors is obtained from symmetry considerations. For example, for the Wannier basis (la, 
Au) of the A = 1 subspace of the diatomic chain, the displacement of the central ion is 
given by one real parameter and each pair of ions equidistant from the center requires an 
additional real parameter to specify the equal displacements (see Fig. ^). The parameters 
of the Wannier basis vector, starting with the innermost coordination shells, are then fit 
to reproduce a set of known normal modes in the band subspace. Since the Wannier basis 
vector is localized, it should be a good approximation to set the displacement parameters 
to zero beyond a finite number of shells. For example, in the diatomic chain, the simplest 
nontrivial truncation is at the first shell. For A = 1, the fit to the eigenvector at A; = yields 
the approximate Wannier basis vectors WiJ with displacements 4= for the central A ion and 
-j= for the first shell of B ions. The frequencies corresponding to the approximate normal 
modes Y,n exp{ikn)wi^ , computed from Eqn. |3.1| and shown in Fig. ^, agree very well with 
the true frquencies if the nonorthonormality of the basis vectors is taken into account. With 
the inclusion of just one more shell of A ions in the approximate Wannier basis vector to 
fit the normalized eigenmode at fc = vr shown in Fig. H, it is nearly indistinguishable from 
the exact Wannier basis vector on the scale of Fig. ^ and unit normalization to within a 
few percent is achieved throughout the Brillouin zone. This construction also illustrates the 
point that the choice of phases for the eigenmodes is important in producing a well-localized 
Wannier basis vector. In this case, if the normalized eigenmode at fc = tt were multiplied 
by —1, the approximate Wannier basis vector would have larger displacements on the first 
A coordination shell than on the central A ion. 

Using the symmetry properties of the Wannier basis vectors, the form of the effective 
Hamiltonian can be written as an invariant expansion of the energy in the Wannier basis 
coordinates of the subspace. We now consider the diatomic chain with effective Hamilto- 
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nian subspace A = 1. Since these coordinates have the same transformation properties as 
individual ionic displacements of the A ions, the effective Hamiltonian parameters can be 
put into direct correspondence with the familiar force constants of the monoatomic chain 
or with tight-binding parameters for a system with one p orbital per unit cell. Truncating 
the expansion at fourth order in onsite terms and at quadratic order and second neighbor 
in intersite terms yields: 

T^effiiUn}) =Yl[ao^An + bo^An\+ Yl «|n-n'|^AnW (3-4) 

" n,n's.t.\n-n'\<2 

The small number of expansion coefficients in this effective Hamiltonian can be determined 
from total energy data for a particular choice of a and 7. Specifically, once the subspace 
Wannier basis vectors {w\n} have been determined, the values of the coordinates {^aw} 
specify a unique ionic configuration. Ti-effii^.An}) is then taken as the potential energy of this 
configuration. Calculation of the coefficients a„ in the invariant expansion proceeds by fitting 
the sum of quadratic terms in Equation |3]^ for coordinates {,^a„ = Re{^k exp{ikn))} , at 
various k, to the quadratic energy of the corresponding ionic configurations, easily obtained 
from the dynamical matrix derived from Equation ^.1| . For example, when ^ = 0.25 the 
quadratic energy at fc = 0, fc = | and k = n can be fit by oq = (13 — \^33)a/l6, ai = —a/4 
and a2 = (—5 + \^33)a/32. Comparison of the quadratic part of the effective Hamiltonian, 
truncated at 02, with the true quadratic energy at other k is used to check the accuracy 
of the truncation of the invariant expansion, which in this case is within a few percent. 
Anharmonic parameters such as 60 are obtained from calculation and fitting of full potential 
energies for finite distortions. 

In the study of structural transitions, focus on a simple description of a relevant mode 
at a particular k-point has occasionally led to a choice of localized basis vectors inconsistent 
with this symmetry criterion. As will be discussed in Section IV, a common example is 
the use of a rotational variable to describe the R25 instability in perovskites. The negative 
consequences of inconsistent symmetry can be simply illustrated in the diatomic chain. The 
displacement pattern of the Xg mode (Fig. H) can be obtained with localized basis vectors 
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not only of the type (16,A„) but also {la,Ag), the latter being inconsistent with the symmetry 
criterion. As shown above, the choice of {lb,Au) is the label of the Wannier basis for the 
A2 band subspace. However, localized basis vectors of the type {la,Ag) lead to zero ionic 
displacements for k = distortions, showing that this "Wannier basis set" is in fact linearly 
dependent. 

The problems encountered when using an inconsistent symmetry choice in the construc- 
tion of the effective Hamiltonian can be understood using Equations |3.2| and |3.3| . Wannier- 
like basis vectors inconsistent with our symmetry criterion can be constructed only if the 
requirement that the inverse transforms (Equation |3.3|) reproduce normalized normal modes 
is lifted and zero normalization allowed. In Equation p.2| , this corresponds to multiplying 
the phase factor exp{i(j)/<^{k)) by a real function fiy{k). For example, Wannier- like functions 
of the type {la,Ag) for the upper branch (A = 2) can then be obtained with 0a(^) = f and 
fh{k) = -. The inverse transforms reproduce the eigenmodes i-e^A, not normalized to one 
except at fc = TT. The energy written in terms of the Fourier transform of the Wannier basis 
coordinates ^a(/c) = J2n^^P{'^kn)C,An is 



E-eAWia(^)P 



k ^ 



Effective Hamiltonian parameters can be fit as above to reproduce the dispersion ^e^ik). 
In particular, at A; = 0, the parameters must be fit to zero energy. To use material-specific 
information at A; = 0, it is necessary to expand the model dispersion to leading order in k"^ 
and to fit the coefficient of -^ to tA{k = 0). The approximate Wannier basis vector con- 
struction can also be carried out as above except at A; = 0, where it is necessary to expand 
the approximate eigenmode to leading order in k and to fit the coefficients to efc=o,A- With 
/a(A; = 0) = 0, required to eliminate the inconsistent condition on the phases at A; = 0, 
the resulting localized basis is linearly dependent and cannot reproduce the normal mode 
at A; = 0. This amounts to eliminating the point k = from the ionic displacement space. 
In the thermodynamic limit, the behavior of the system should not be affected, since in 
the partition function configuration integral, this is a set of measure zero. However, in a 
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finite-size simulation, the integral is replaced by a sum, and the incorrect description of 
the system arising from the linear dependence leads to an 0{j^) error. If transformation 
properties consistent with our symmetry criterion are chosen and the normalization require- 
ment imposed, these finite-size simulation errors and the unnecessary complications in the 
first-principles fitting procedure are avoided. 

IV. EFFECTIVE HAMILTONIANS FOR CUBIC SYSTEMS 

In this section, we generalize the construction of the effective Hamiltonian to arbitrary 
structures and apply it to phonon-related structural transitions in the perovskite, rocksalt, 
fluorite and A15 structures. We will use as our primary illustration the example of the 
perovskite structure. Among the many compounds in this structure, a rich variety of phonon- 
related structural transitions are found.! These can be analyzed in a unified systematic 
framework using the present approach. 

The cubic perovskite structure, which will be taken as the high-symmetry reference struc- 
ture, has five atoms per unit cell and space group PmSm. The first step in the effective 
Hamiltonian construction is a symmetry analysis to identify the band subspaces (i.e. mini- 
mal invariant subspaces spanned by one or more entire branches of normal modes). Once the 
space group irrep labels for the ionic-displacement space have been determined, a complete 
list of the band subspaces and the corresponding Wannier basis labels can be obtained by 
the following general procedure, valid for arbitrary structures. First, all possible symmetry 
types of band subspaces are constructed using the compatibility relations between the irrep 
labels at the high-symmetry fc-points. The corresponding Wannier basis labels are found by 
considering in turn eachWyckoff symbol with each irreducible representation of the corre- 
sponding site symmetry group such that the product of the number of the Wyckoff positions 
with the dimension of the site symmetry group irrep is less than or equal to the maximum 
multiplicity of the band subspaces. A particular combination will label a Wannier basis of a 
band subspace if the space group irrep labels produced at the high-symmetry fc-points match 
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those of the band subspace.El Apphcation of this procedure to the perovskite structure yields 
four band subspaces, all labelled by occupied Wyckoff positions with vector representations, 
corresponding to atomic displacements. These are listed in Table I. 

As in the case of the diatomic chain, a parametrization of the displacement pattern 
of the Wannier basis vectors can be obtained from symmetry considerations. For each 
band subspace, the procedure is (1) choose one of the Wyckoff positions specified by the 
label; (2) identify the symmetric coordination shells surrounding that site; and (3) compute 
the independent displacement patterns of each shell that transform according to the given 
irreducible representation of the site symmetry group. For each of the four band subspaces 
of the perovskite structure, the displacement patterns for the innermost coordination shells 
are shown in Fig. p. The relative amplitudes of each independent displacement pattern for 
each shell are free parameters determined by the requirement that a particular branch or 
branches of normal modes be reproduced, and thus will be different for individual materials. 

As we have seen in Section II, construction of an effective Hamiltonian involves the choice 
of an invariant subspace as the direct sum of one or more band subspaces. The effective 
Hamiltonian subspace should include all modes relevant to the particular transition being 
modelled. The smallest subspace satisfying this condition, and therefore the maximum 
possible reduction in the degrees of freedom, can be identified by listing the symmetry 
labels of the relevant modes and choosing a direct sum of the band subspaces of the ionic- 
displacement space which produces the desired labels. There are cases in which there is 
more than one minimal choice which satisfies the conditions. If the choices differ in the other 
symmetry labels, the choice is made for an individual material by determining the symmetry 
labels of the lowest-energy modes at other high-symmetry A;-points. If all symmetry labels are 
the same, then the Berry phase labelc3£3 determines the choice, since well-localized Wannier 
basis vectors are only obtained by transforming eigenvector branches which are smooth in 
k. The general analysis of Berry phases for composite energy bands is quite complicated.Ea 
For present purposes, we observe that if only one of the possible choices is produced by an 
occupied Wyckoff position combined with a vector representation, that choice will smoothly 
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reproduce the physical eigenvector branch. Examples of systems where more than one choice 
is of this kind can be constructed, the simplest being a crystal of space group F222 (^ 22)EZl 
with both the 2a and 2b Wyckoff positions occupied. In such cases, it is advisable to include 
both branches in the effective Hamiltonian subspace, with two Wannier basis vectors, one 
centered on each Wyckoff position. Only if the "relevant phonon" is heavily dominated by 
displacements of only one of the two atom types should construction of a subspace with only 
one Wannier basis vector (centered on the corresponding Wyckoff position) be considered. 
In the crystal structures analyzed in this section, this difficulty does not arise. 

With the selection of the effective Hamiltonian subspace, the invariant expansion of 
the lattice energy in Wannier basis coordinates proceeds straightforwardly, facilitated by 
the symmetry properties of the Wannier basis vectors. Different structural transitions ob- 
served in perovskite structure compounds require different treatment, starting at this point. 
In the following, we will consider four common perovskite transitions (i) the ferroelectric 
transition in BaTiO^ and PbTiO^, with an unstable Fis mode; (ii) the oxygen-octahedron 
rotational transition in SrTiO^, with an unstable i?25 niode; (iii) the Jahn- Teller instabil- 
ity in Lai-x{Ca, Sr, Ba)xMn03, with an unstable R2 mode; and (iv) the antiferroelectric 
transition in PbZrOs, involving several unstable modes. 

For the Fis ferroelectric transition found in BaTiO^ and PbTiO^, all four subspace 
choices A, B, O^^x and O^^y are consistent with the symmetry criterion. The simplest form 
of the effective Hamiltonian is obtained by choosing either the 3A^-dimensional subspace A, 
resulting in the invariant energy expansion for PhTiO^ of Ref. ^, or the 3A^-dimensional 
subspace B for BaTiO^ as in Ref. p. In these two cases, the choice of A versus B is 
determined by the symmetry labels of the lowest zone-boundary modes. In these systems, 
coupling to strain is known to be important.H'Ea'Ea Homogeneous strain is easily included 
by introducing the six-component strain variable 6^/3 (a,/5 = x,y,z) and its lowest order 
coupling to the ferroelectric mode coordinates.cra To include inhomogeneous strain as well, 
the subspace can be enlarged by adding a second band subspace containing F15 to describe 
the acoustic modes. The most convenient choice is the cation subspace not already used 
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for the ferroelectric mode. In addition to the space group symmetries, the inhomogeneous 
strain Hamiltonian has global translation and rotation invariances which can most easily be 
satisfied by following the construction method of Keating.E3 If the expansion is truncated 
at three independent parameters, corresponding to the three elastic constants, the strain 
Hamiltonian used in Ref. ^ is obtained. Additional expansion parameters could be included 
to improve the description of the acoustic modes away from the zone center. With Wannier 
basis vectors that exactly reproduce the mode eigenvectors, there is no quadratic coupling 
between the two subspaces A and B. The lowest order anharmonic coupling that does not 
vanish in the limit k —>■ is the coupling term used in Ref. p. 

The i?25 oxygen-octahedron rotational transition occurs in SrTiO^ at 105 K and in a 
number of other systems.Ell This unstable mode is built up in Refs. ^ and from a Wannier- 
like "rotational variable" transforming according to TipE3 at each site of the simple cubic 
lattice. Symmetry analysis shows that this basis set is linearly dependent, generating sym- 
metry labels r25, Ag and X2 which do not appear in the symmetry decomposition of the 
ionic displacement space. While this difficulty can be formally resolved by requiring the 
energy to be zero for configurations with /c in F — A — X, the vanishing of the normalization 
at lines in the Brillouin zone leads to 0{^) errors in finite-size simulations, as discussed 
in the previous section for the diatomic chain. Also, incorporation into the effective Hamil- 
tonian of first-principles results at these lines requires extra manipulations. Finally, use of 
this basis in quantum-mechanical or dynamical studies, in which the kinetic energy must 
be treated explicitly, would complicate the analysis by introducing an artificially vanishing 
effective mass in the neighborhood of these lines. 

In the present approach, to describe the i?25 oxygen-octahedron rotational transition 
in SrTiO^, the band subspace is the 6A^ dimensional space O^^y The resulting effective 
Hamiltonian describes a system of two-component vectors on the faces of the unit cells of 
the simple cubic lattice (Fig. |3(a)). This form is slightly more complicated than the models 
usually studied in the statistical mechanics literature, since the "spin" sites do not form 
a Bravais lattice. However, symmetry still ensures that the invariant energy expansion is 



fairly simple. In fact, the quadratic part has the form of a submatrix of the familiar force 
constant matrix for the perovskite structure. Furthermore, this form presents no problems 
for numerical analysis of the model through Monte Carlo simulation. Even the fact that the 
band subspace is six dimensional can be turned to advantage. The extra Fis mode can be 
used to describe inhomogeneous strain coupling or the ferroelectric mode, which is observed 
to be nearly unstable.E^ 

To describe the i?2 Jahn- Teller instability in Lai^^iCa, Sr, Ba)xMnO^ the invariant 
subspace for effective Hamiltonian construction must include the 3A^ dimensional band sub- 
space labelled Ox,x in Table I. The resulting effective Hamiltonian describes a system of 
one-component vectors on the edges of the unit cells of a simple cubic lattice. As in the case 
of SrTiOs, the invariant energy expansion turns out to be simple and quite tractable for 
analysis through Monte Carlo simulation. The symmetry of the model can also be readily 
exploited in writing the form of the coupling to electronic states.la 

For antiferroelectric PbZrO^, there are four relevant modes: an R25 oxygen-octahedron 
rotation, a F15 polar mode, a (||0)(^)S3 mode, and a (110) (^)M5 mode.E3 All these modes 
are included in the 6A^-dimensional band subspace Ox,y (Fig. 0(a)), so that the resulting 
effective Hamiltonian is surprisingly simple for such a complex transition. To examine the 
effects of coupling to inhomogeneous strain, the effective Hamiltonian subspace would have 
to be extended to include an additional 3A^-dimensional subspace [A or B). 

The analysis for structural transitions occurring in other structures proceeds similarly. 
The rocksalt-rhombohedral transition in GeTe has been previously studied with a first- 
principles model Hamiltonian approach.Q The high-symmetry reference structure is the rock- 
salt structure, with an fee Bravais lattice with conventional lattice constant oq and space 
group FmSm. The symmetry analysis of the band subspaces is given in Table II. The soft 
mode of the transition to the rhombohedral structure is F15, so both Ge and Te are possible 
choices for the effective Hamiltonian subspace. As in the perovskite structure, the choice is 
determined by the symmetry label of the lowest energy zone-boundary mode. The effective 
Hamiltonian constructed with the Te subspace is a system of three-component spins on an 
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fee lattice of lattice constant oq. In Ref. ^, coupling to homogeneous strain was included and 
good agreement obtained with the experimental transition behavior. The rocksalt structure 
is so simple that if inhomogeneous strain were included in the model (by adding the Ge sub- 
space) the resulting "effective Hamiltonian" subspace would be the full ionic displacement 
space, and the effective Hamiltonian would be identical to Hiat- 

The cubic-tetragonal transition in Zr02 has inspired several recent first-principles band- 
structure and total-energy studies.E3 In this system, the high-symmetry reference structure 
is the fluorite structure, with an fee Bravais lattice with conventional lattice constant Oq and 
space group Fm3m. The symmetry analysis of the band subspaces is given in Table III. The 
soft mode of the transition to the tetragonal structure is X2, so the effective Hamiltonian 
subspace must be the 6A^-dimensional space O. The resulting effective Hamiltonian is a 
system of three-component spins on a simple cubic lattice of lattice constant y- However, 
the symmetry for the invariant energy expansion is not Oh, but T^, imposed by the presence 
of the fee zirconium sublattice. 

Structural transitions in A15 compounds such as V^Si have long been of interest because 
of the interplay with "high-temperature" superconductivity in these systems. E3 The A15 
structure has two formula units (8 atoms) per unit cell, with a simple cubic Bravais lattice 
and a nonsymmorphic space group Pm3n. The symmetry analysis of the band subspaces is 
given in Table IV. The observed low-temperature structure is generated by the ri2 mode. 
The effective Hamiltonian subspace is then the 6iV-dimensional space Ay. The choice of 
y4|| over 6d,u\i^, which also produces the right symmetry label, is made on the basis of the 
Berry phase argument discussed earlier in this section. The resulting model is a system of 
one-component vectors on three interpenetrating sets of chains running along the cartesian 
directions. The Fis mode in this subspace can be used to include the effects of coupling to 
inhomogeneous strain. 
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V. PRACTICAL ASPECTS OF FIRST-PRINCIPLES IMPLEMENTATION 

The information available from first-principles methods is (i) total energy of supercells; 
and (ii) force constant matrices at various k calculated using density functional theory linear 
response.Ejila In this section, we present a practical approach for using this information to 
construct the effective Hamiltonian for a particular system, illustrating the procedure with 
the application to the ferroelectric transition in PhTiO^,^ 

The analysis begins with the identification of the high-symmetry reference structure and 
relevant phonon, and the first-principles calculation of the force constant matrices (either by 
linear response or by frozen phonon supercell calculations) at the high-symmetry /c-points. 
The first step is to identify the effective Hamiltonian subspace Aq. For some structural 
transitions, Aq will be a single band subspace uniquely specified by the relevant phonon. 
For others, the relevant phonon allows a choice of band subspaces. In the example of the 
PhTiO^ transition, the Fis phonon allows a choice of band subspaces A, B, Ox,x and O^^y 
The choice should be made so as to include the lowest energy eigenmodes at other high- 
symmetry A;-points. In PbTiO^, the low energy modes X'^, Mg, M^, and R15 determine the 
efffective Hamiltonian subspace to be the band subspace A (see Table |). 

The next step is the construction of the Wannier basis vectors for the subspace Aq. 
The strategy is to fit independent displacement parameters for the innermost coordination 
shells to the normalized eigenvectors at only the high-symmetry /c-points (with adjustment 
of the eigenvector phases to obtain optimally rapid decay). Since the Wannier basis vector 
ionic displacements generally decay rather rapidly with distance from the central ion, this 
is sufficient to determine an excellent approximation to the Wannier basis vector. In the 
example of the diatomic chain, discussed in Section HI, it was sufficient to fit only the k = 
and k = n modes. The application to ionic insulating crystals, including PbTiO^, requires 
additional comment. If Aq contains a degenerate polar A; = mode and there are polar k = 
modes not included in Aq, (as, for example, in the perovskite structured), the electric field at 
k -^ mixes the LO modes differently from the TO modes, and a single Wannier basis vector 
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cannot faithfully reproduce both the transverse and longitudinal branches. Therefore, the 
normal modes of the corresponding longitudinal branch should not be used in constructing 
the approximate Wannier basis vector. The correct approach is to use only the low energy 
(transverse) branches to determine the approximate Wannier basis vector. For PbTiO^, 
we use displacements of the central, first and second neighbor Pb ions and the nearest Ti 
and O shells, for a total of nine independent parameters (six of which are included in Fig. 
P) to reproduce the normalized lowest-frequency normal modes Fis, i?i5 and Mg and the 
normalization of the M5 mode exactly. The amplitudes of the displacements are found 
to decay rapidly, as expected. Moreover, the X'^ and M5 normal modes predicted from 
the parameters obtained from the fit are in excellent agreement with the corresponding 
normalized normal modes calculated from first principles.S 

The form of the effective Hamiltonian is obtained by Taylor expanding the potential 
energy in symmetry invariant combinations of the C,i and approximating to obtain an ex- 
pression with a relatively small number of effective Hamiltonian parameters which can be 
determined from the available information. For the quadratic part of the effective Hamilto- 
nian, this will mean that intercell interactions of the most general form are included only up 
to a finite range. For PbTiO^, we include general intercell interactions up to third neighbor 
(eight independent parameters). The long-range dipolar interactions required for the correct 
description of polar modes can be included in a parametrized form, requiring only knowledge 
of the Born effective charges Z* and the dielectric constant Coo, which can be obtained from 
separate linear response calculations or from energies of distortions with q close to zero. 
Similarly, a small number of higher order terms must be selected so as to ensure stability 
and a ground state configuration of the correct type. In PbTiO^, only onsite anharmonic 
terms are needed. Lastly, in a cubic system such as PbTiO^, inclusion of the coupling to 
homogeneous strain introduces only a few additional parameters. 

The calculation of the parameters in the effective Hamiltonian exploits the fact that a 
given set of Wannier basis coordinates uniquely determines the positions of the ions. The 
quadratic parameters are obtained by fitting to the frequencies of a set of modes in the ef- 

22 



fective Hamiltonian subspace computed from first-principles force-constant matrices at the 
high-symmetry points A;-points and at A;-points along selected symmetry lines. In PbTiO^, we 
compute Z* and eoo from DFT perturbation theory, and fit the eight short-range interaction 
parameters to the following effective Hamiltonian eigenmodes: {Ti^, Rn^, M2, M^, X2, X'^} 
and the A3 modes at ^^(lll) and ^(111). The truncations in the Taylor expansion of 
the quadratic effective Hamiltonian can be checked by computing force constant matrices 
at additional k. With linear response, this overdetermination of parameters is relatively 
easy, while it can be quite impractical in supercell calculations. In the same way, anhar- 
monic terms can be obtained from fitting the higher-order dependence of first-principles 
total energies on the amplitude of the Wannier basis coordinates. Lastly, the homogeneous 
strain terms are obtained by computing the dependence of total energies of configurations 
on variations in the lattice parameters. 

The approximate Wannier basis vector approach can be applied successfully even if 
the branches of interest are not fully isolated in energy from other branches, which are 
recognizable by distinct symmetry labels. Even in this case, a simple model, with the 
effective Hamiltonian acting in the Aq subspace, should still give accurate results for finite- 
temperature properties if the A 7^ Aq branches cross only the higher-energy portions of the 
Aq branches, preferably at A;-points far from the k of the relevant phonon. In that case, the 
approximate Wannier basis vector is fit using only the normalized eigenvectors at A;-points 
where the Aq branches can be distinguished from the A 7^ Aq branches by their symmetry 
labels. This set of A;-points will typically be the same as the set of high-symmetry A;-points 
prescribed in the general procedure for approximate Wannier basis vector construction. 
Similarly, at other /c-points, energies of the approximate eigenvectors built up from these 
approximate Wannier basis vectors are computed from force constant matrix results and used 
to determine effective Hamiltonian parameters. Therefore, application of the approximate 
Wannier basis vector approach essentially projects out the A 7^ Aq character at the low- 
symmetry A;-points, which amounts to ignoring the crossing. 

If a substantial amount of crossing occurs into the lower energy range, if all the relevant 
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modes are not contained in a single band subspace, or if anharmonic coupling to modes 
outside the relevant band subspace is important, it will be necessary to enlarge the effective 
Hamiltonian subspace to include the additional branches by working with a direct sum 
Ai © A2 © ... © A^. For each band subspace, there is an independent Wannier basis vector 
to be determined by fitting eigenvector information. If the approximate Wannier basis 
vectors exactly reproduce all branches at all calculated /c-points, then in the invariant energy 
expansion there are no cross terms between coordinates from different band subspaces. If 
this is not the case (for example, for the polar modes discussed above) cross terms must be 
included in the invariant energy expansion and explicitly determined from the force constant 
matrix information. 

If the subspace containing the relevant modes contains no polar modes and is well sepa- 
rated in energy from the higher-energy eigenmodes, there is an alternative approach which 
bypasses the explicit construction of the Wannier basis vectors altogether by applying the 
philosophy of tight-binding parameterization of electronic energy bands.cH In that case, we 
assume the existence of underlying Wannier basis vectors and fit directly to the eigenvalues 
to obtain effective Hamiltonian parameters. This approach was discussed in Section III for 
the diatomic chain example, which satisfies the necessary conditions. In applications to real 
materials, this ideally simple situation is not expected to be very common, and therefore 
the approximate Wannier basis vector construction is recommended for routine use. 

VI. CONCLUSIONS 

A systematic procedure has been presented for the construction of localized Wannier 
basis vectors and corresponding effective Hamiltonians, applicable to arbitrary structures 
and phonon-related transitions. The application to several different structures has been 
described. For the perovskite structure, a new approach to the oxygen-octahedron tran- 
sition in SrTiOs is suggested, which avoids the subtle difficulties arising from the use of 
a "rotational variable." Improvements can be systematically made by including additional 
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band subspaces in the effective Hamiltonian subspace, if indicated by the presence of sub- 
stantial crossing of eigenmode branches or strong anharmonic couphng to inhomogeneous 
strain or other modes. The effective Hamihonian parameters for a particular material can 
be obtained using any method that can reliably calculate force constants and distortion 
energies. Once the force-constant matrices have been calculated, the construction of an 
effective Hamiltonian is straightforward and requires very little additional computational 
effort. An efficient method, such as DFT linear-response, expedites the necessary overdeter- 
mination of effective Hamiltonian parameters. Further applications of this approach could 
include the construction of realistic quantum mechanical models for the study of quantum 
paraelectrics and ferroelectrics,E3 as well as classical dynamical models for the study of ferro- 
electric phenomena such as switching. In addition, these effective Hamiltonians can be used 
for realistic first-principles studies of thermal expansion and other temperature dependent 
lattice properties,c3 as well as for structural transitions. More generally, the construction of 
a Wannier basis may prove advantageous for the treatment of localized perturbations. 
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TABLES 
TABLE L Band subspaces of the ionic displacement space of the perovskite structure. Sub- 
spaces are labelled by the corresponding Wyckoff symbol and site symmetry group irrep. In cases 
where the subspace is generated by a set of occupied Wyckoff positions and a vector representation, 
the atom type is used as the label. For each subspace, the symmetry labels at the high-symmetry 
fe-points are as in Ref. 23, with conventions taken from Ref. E3 and origin at A. 
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TABLE IL Band subspaces of the ionic displacement space of the GeTe rocksalt structure, as 
in Table L Symmetry label conventions are taken from Ref. ^6], with the origin at Te. 

r L X W 

Ge Lis Li,L3 X'^,X'^ W^W^ 

Te ri5 L'2,L'3 X'^,X'5 W^,W3 

TABLE in. Band subspaces of the ionic displacement space of the Zr02 fluorite structure, as 



in Table I. 
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TABLE IV. Band subspaces of the ionic displacement space of the A15 structure, as in Table I. 
^11 denotes an A ion displacement along the chain direction and A± denotes an A ion displacement 
perpendicular to the chain direction. For the 66 sites, the irreducible representations of the site 
symmetry group are labelled by the component of a polar vector (n) or an axial vector R parallel to 
the chain (||) or perpendicular and pointing along the direction to the nearest B atom (_L, S). For 
the 6d sites, the polar vector components are parallel or perpendicular to the "chain" directions 



defined by these sites. Symmetry labels at F are taken from Ref. 36 and at zone boundary points 
from Ref. ^ with the origin at B. 
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FIGURES 
FIG. 1. The harmonic potential energy of the diatomic chain model is shown by two types of 

springs. One type, with spring constant a, connects nearest neighbor A ions (open circles) and B 

ions (solid circles). The other, with spring constant 7, connects nearest neighbor B ions. 

FIG. 2. Eigenvalues are shown as a function of k for the two branches of normal modes of the 
diatomic chain with — = 0.25. At A; = and k = n, each branch is labelled by the appropriate 
irreducible representation of the space group and the corresponding displacement pattern is shown 
in the adjacent sketch. A ions are shown as open circles and B ions as solid circles. 

FIG. 3. The harmonic potential energy of the triatomic chain model is shown by two types of 
springs. One type, with spring constant a, connects nearest neighbor A ions (open circles) and B 
ions (solid circles). The other, with spring constant 7, connects nearest neighbor B ions. 

FIG. 4. The displacement pattern of the Wannier function w\n, corresponding to the lower 
branch in Fig. ^, is shown by the arrows. The pattern is odd around the central A ion. The arrows 
for the outermost A and B ions are smaller than the ion symbols. 

FIG. 5. Dispersion '^A., for the lower branch computed using the approximate Wannier func- 
tion fit to A; = (lower solid curve) compared with the true dispersion (dashed curve). The upper 
solid curve shows the overlap matrix element S" as a function of k. For comparison, the upper 
dashed line is the unit normalization S = 1. 

FIG. 6. Displacement patterns for the innermost coordination shells of representative Wannier 
basis vectors in each of the four band subspaces of the perovskite structure. The open squares 
represent the A ions, the shaded squares the B ions, and the open and shaded circles represent 
inequivalent oxygen ions. The four lines in the figure correspond to the subspaces A, B, Ox,x and 
Ox,y In the first structure in each line, the displacement of each set of inequivalent ions can be 
specified independently. 
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FIG. 7. (a) The model system for the i?25 oxygen-octahedron rotational transition in SrTiOz 
consists of two-component vectors on the faces of the unit cell of a simple cubic lattice; (b) the model 
system for the R'2 Jahn- Teller instability in Lai-x{Ca,Sr^Ba)xMnO^ consists of one-component 
vectors on the edges of the unit cells of a simple cubic lattice. 
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